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Abstract 

We investigate numerically the long-time behavior of balanced 
Alfven wave turbulence forced at intermediate scales. Whereas the 
usual constant-flux solution is found at the smallest scales, two new 
scalings are obtained at the forcing scales and at the largest scales of 
the system. In the latter case we show, in particular, that the spec- 
trum evolves first to a state determined by Loitsyansky invariant and 
later a state close to the thermodynamic equipartition solution pre- 
dicted by wave turbulence. The astrophysical implications for galactic 
magnetic field generation are discussed. 

1 Introduction 

Turbulence flows are ubiquitous in astrophysical environments from the solar 
wind (Goldstein et al., 1999; Galtier, 2006), to interstellar (Elergreen et al., 
2004; Scalo et al., 2004), galactic and even intergalactic media (Govoni et al., 
2006). At the larger scales, signatures of astrophysical turbulence are found, 
in particular, in the magnetic fleld measurements whose origin remains one 
of the major challenging problems (Pouquet, 1993; Widrow, 2002; Branden- 
burg, 2005). In this paper, we emphasize a new mechanism for generating 
large-scale magnetic fleld. It is based on the resonant interactions of shear- 
Alfven waves in a turbulent medium permeated by a strong external magnetic 
fleld and influenced by an external forcing at intermediate scales. In this sit- 
uation, both large-scale kinetic and magnetic energies may be produced by 
mainly non local interactions. This scenario, although very simple, may be 
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relevant for describing re-generation (i.e. maintenance) of large-scale galactic 
fields, e.g. in our galaxy where energy is injected at intermediate scales by 
stellar winds and supernovae explosions on scales 10 — lOOpc (Ferriere et al., 
2004). 

2 Alfvenic turbulence 

The wave kinetic equations of incompressible MHD were derived rigorously 
by Galtier et al. (2000,2002). Let us consider the simplest case when helicities 
are taken to be zero, when pseudo-Alfven waves are discarded and when 
turbulence is axially symmetric and balanced (equal spectra for Alfven waves 
CO- and counter-propagating with respect to the external magnetic field). In 
this case for the shear-Alfven wave energy spectrum we have E{k±,k\\) = 
f{k\\)E±{k_i), where /(A;||) is a function fixed by the initial conditions (i.e. 
there is no energy transfer in the parallel direction), and the transverse part 
obeys the following kinetic equation, (Galtier et al., 2000) 

dtE±(k±) + uklEAkA = ^ [ cos^ (f)sme —EAk±) 

[k±E^{L^) - L^E^ik^)] dK^dL^. (1) 

Here feo the background magnetic field (which will be taken equal to one), 
(p the angle between k_|_ and L_|_, and 9 the angle between k_|_ and k±. This 
integro-differential equation has been computed in (Galtier et al., 2000) but 
only for short times in order to analysis, in particular, the formation of the 
Kolmogorv-Zakharov solution. In this paper we investigate the long-time 
behavior of Alfven wave turbulence when of a forcing term is applied. In 
principle, this regime may be investigated for any value (7^ 1) of the magnetic 
Prandtl number; in this case the viscosity v has to be seen as a combination of 
the kinematic viscosity and the magnetic diffusivity. Note that in the alfvenic 
turbulence regime, the kinetic and magnetic energies are equal, therefore all 
results may be interpreted in terms of magnetic energy fluctuations. 

2.1 Short-time evolution 

We perform numerical simulations of the wave kinetic equation ([1]) of MHD 
turbulence in which a forcing term is introduced at intermediate scales. A 
non- uniform grid in Fourier space is used (Galtier et al., 2000) with k^{i) = 
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Figure 1: Evolution of the energy spectrum for times = (solid), ti = 
5 X 10""^ (long dash), t2 = 2ti (small dash), ^3 = 3ti (dot-dash), ^4 = Ati 
(three dots-dash) and ^5 = 5ti (bold solid). For comparison, three slopes are 
given with power law indices -|-3, —1.4 and —2 respectively. 



Sk 2'/^ {i e [0,156]; 5k = 1/8; F = 8). This method allows us to achieve 
very high Reynolds numbers (> 10^ with u = t] = 5 x 10^^). A weak forcing 
is applied to the flow at wavenumbers k± E [2^^^^, 2"*^^/^] with a constant 
(flat) spectrum. 

In Figure [1] the evolution of the shear-Alfven wave energy spectrum is 
displayed for short times. The forcing scales determine the initial evolution 
as can be seen at intermediate wavenumbers. The first robust property to be 
observed is the generation of small-scales through a direct cascade and the 
formation of the well-know kj^ solution predicted by wave turbulence theory 
(Galtier et al., 2000). At the same time we see that the largest scales of the 
MHD flow are affected by the forcing and produce a power law slightly steeper 
than the reference slope (-1-3). At intermediate scales, no clear tendency is 
observed yet, only a reduction of the flat spectrum seems to occur. 

In Figure [H we show the temporal evolution of the transfer function. 



3 



T{kj_), and the energy flux, P{k±), for tlie same times as in Figured! We re- 
mind tlie general relation in the inertial range (where forcing and dissipation 
are negligible) 

a,Bx^r^-|l. (2) 

The constant flux solution of wave turbulence is clearly observed at small- 
scales and corresponds to a direct cascade (positive flux). As expected, a 
negative transfer is found at the forcing scales but a strong asymmetric proflle 
is obtained initially which reduces with time. This reduction of asymmetry 
appears when the fcj^ solution is well established (see Figure [T]). At the 
largest scales a positive transfer with a negative (non constant) flux are 
found which means that energy is accumulated probably mostly by non-local 
interactions. Note that this is not an inverse cascade behaviour because there 
is no extra (in addition to the energy) positive invariant in our problem that 
could lead to a dual cascade process. 



2.2 Intermediate-time evolution 

In Figure [3], we show a next set of times that we call intermediate-time. 
Whereas the small-scale power-law spectrum appears to be stable with no 
apparent change, at the largest scales a scaling is found. This solution 
could be related with the invariance of the Loitsyansky integral (Davidson, 
2004; Bigot et al., 2007) which leads to scaling for the modal spectrum in 
at low wavenumbers with s = D + 1 [D being the space dimension). This 
analysis is often advocated to justify the scaling taken initially in (direct) 
numerical simulations. For example, it is shown in a recent work (Bigot et 
al., 2007) that the k^ scaling is well conserved in freely decaying wave MHD 
turbulence. At the intermediate scales, we see clearly a new scaling close to 
fcj^'^. It takes about one order of magnitude longer to build this solution than 
to form the small-scale cascade spectrum. We will see that this solution may 
be explained rigorously by an analysis based on the wave kinetic equation. 

Figure H] gives the transfer function and flux for the same times as in 
Figure [3l The constant flux at small-scale is still observed with an increase of 
its value in time; note that this increase seems to saturate. At wavenumbers 
around one, we see that the transfer function exhibits a change of sign going 
from a positive value, for wavenumbers smaller than one, to a negative value, 
for wavenumbers larger than one, which traduces an energy transfer towards 
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Figure 2: Temporal evolution of the transfer function T{k±) and the flux 
P{k±) at the same times as in Figur§[T] (the same notation are also used). 
Top: general behavior at all scales of T and P; note the formation of the con- 
stant flux solution at small-scales. Bottom: zoom at large-scales; a negative 
transfer is clearly seen at intermediate (forcing) scales whereas a negative 
flux is found at the largest scales. 
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Figure 3: Evolution of the energy spectrum for times = 0.75 x 10~^ (solid), 
t7 = 1.75 X 10-2 (long dash), tg = 2.75 x lO'^ (small dash), tg = 3.75 x lO'^ 
(dot-dash), tio = 4.75 x 10"^ (three dots-dash) and tu = 0.1 (bold sohd). 
The same slopes as in Figure [T] are given. 
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Figure 4: Temporal evolution of the transfer function T{k±) and the flux 
P{k±) at the same times as in Figure'El Top: general behavior at all scales. 
Bottom: zoom at large-scale. 
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Figure 5: Long-time evolution of the energy spectrum at ti2 = 0.15 (solid), 
ti3 = 0.35 (long dash), = 0.55 (small dash), ti5 = 0.75 (dot-dash), tig = 
0.95 (three dots-dash) and tu = 1.45 (bold solid). 

larger scales. In the meantime, we note at the largest scales a decreasing of 
the flux which is still negative. 

2.3 Asymptotically long-time evolution 

The long-time behavior of Alfven wave turbulence has been investigated as 
well. Figure [5] shows the energy spectrum evolution for times up to t = 1.45. 
Only the large-scale flow evolves at this stage, namely from the scaling 
to a flatter power law. 

Figure [6] gives more interesting information. The tendency to flux satura- 
tion observed at the smallest scales is confirmed and the transfer function at 
intermediates scale does not change anymore. As expected, only the largest 
scales are evolving with a transient decrease of the absolute value of the flux 
(which is negative at the large scales) accompanied by a plateau formation 
and a decrease of the transfer function. 

The last set of graphs displays the final behavior (up to time t = 24.5) of 
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Figure 6: Temporal evolutions of the transfer function T{k±) and the flux 
P{k±) at the same times as in Figure^E 
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Figure 7: Asymptotically long-time evolution of energy spectrum at times 
tis = 2.5 (solid), ti9 = 4.5 (long dash), ^20 = 6.5 (small dash), ^21 = 8.5 
(dot-dash), t22 = 10.5 (three dots-dash) and t23 = 24.5 (bold sohd). 

Alfven wave turbulence. Figure [7] reveals the formation of a new power law 
close to /c^^ at the largest scales whereas the other scales are still steady. 

In Figure M the flux weakening and the plateau formation is continuing 
at the largest scales whereas the transfer function is still decreasing. 

3 Discussion 

3.1 Nonlinear dynamics at intermediate (forcing) scales 

The scaling behavior found at intermediate scales seems to be strongly de- 
pendent on the type of forcing applied to the MHD flow. This is confirmed 
with a second simulation (not shown) where a non fiat forcing is applied. The 
nonlinear dynamics at forcing scales may be investigated through a rigorous 
analysis based on the wave kinetic equations. Let us start from equation ([1]) 
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Figure 8: Temporal evolutions of the transfer function T{k±) and the flux 
P{k±) at the same times as in FigurdUl 



and assume a power law scaling solution 



E^{k^) = Ckl. (3) 
Then let us assume the following general forcing 

fik±) = Cjkl (4) 

Therefore, we have 

dtE^iki.) + uklE^ik^) = C.I. + fik^), (5) 

where C.I. is the collision integral. We focus our analysis on the forcing scales 
and neglect the dissipative term. By introducing the notation k± = k±x and 
L± = k±y, we obtain the steady solution 

C^k^^+^I{n) - Cfki, (6) 

where the collision integral writes 



Kn) 



,2,1 2 



+ 1 — X 



2 _i_ 1 2 \ 

^ ~^ ^ ~ y \ n-l_/_n-l 



Since I{n) is not singular, we directly obtain the scaling relation 

n=^- (8) 

For a constant (flat spectrum) forcing, we get n = —3/2. The solution found 
is indeed very close to this power law solution. The slight difference observed 
may be attributed to the non exact stationarity of the solution (see Figure [9]). 
Note that we have tested the theoretical prediction against another forcing 
with ^ = 1. We find a steady energy spectrum in kj^'^ at intermediate scales 
which is close to the prediction in kj^ whereas the large-scale spectrum seems 
not to be affected and behaves similarly as in the simulation presented in this 
paper. 



3.2 Solution at the largest scales 

The asymptotically long-time evolution is characterized at the largest scales 
by a decreasing and fiat in k energy flux (see Figure [HD. Contrary to what 
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happens at the smallest scales, the flux has a very small value (around —5 
compare to 4x 10^ at small scales) whereas the transfer function tends to zero. 
Therefore, the flnal state of Alfven wave turbulence seems to be the zero-flux 
thermodynamic solution which corresponds to equipartition of energy in the 
/c-space. In this case the theoretical solution derived in (Galtier et al., 2000) 
gives 

E ~ k^. (9) 

We see that the scaling found numerically is somewhat flatter, i.e. k^^. The 
difference may be attributed, in particular, to the limited time of computation 
which does not allow to reach exactly the zero-flux solution. The problem 
is that the evolution slows down drammatically as one moves closer to the 
lowest wavenumbers, and it takes an incredibly long computational time to 
reach the flnal state. 

3.3 Global behavior 

To complete our analysis, the temporal evolution of total shear-Alfven wave 
energy is displayed in FigureOfor times up to t = 2.5. Different time scalings 
are clearly visible: a flrst power law in t is found and a second in about t^/^ 
is estimated. The flrst temporal power law is linked to the formation of 
the Kolmogorov-Zakharov solution during which a transient solution with 

'V In 

spectrum k_^ is produced (Galtier et al., 2000). The second temporal 
power law corresponds to the formation of the steady spectrum solution at 
intermediate scales discussed above. Finally, after time t > 0.2 a plateau 
is reached with negligible energy variation in time (at least in logarithmic 
coordinate) . 

4 Conclusion 

We have seen that the large-scale magnetic fleld is produced on a time-scale 
which is about 10^ larger than the time-scale needed to establish the small- 
scale energy cascade solution. We believe that the flnal distribution in the 
large-scale part should correspond to the thermodynamic energy equiparti- 
tion k]_ even though it would take an extremely long computing time for 
formation of this spectrum which we were not able to achieve. Instead, we 
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Figure 9: Temporal evolution of shear- Alfven wave energy for times up to 
^18- For comparison, two slopes are given with power law indices +1 and 
+1/3 respectively. The vertical lines correspond to times ti to tn. 
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observed formation of a shallower long-term scaling at the large scales, kj^, 
which we believe to be transient. 

Because of the energy cascade to the smallest dissipative (resistive and 
viscous) scales, the small-scale part of the spectrum is important for un- 
derstanding the total energy dissipation rate. On the other hand, it is this 
large-scale part of the spectrum that contains most of the wave (magnetic 
and kinetic) energy itself. One can view it as a sort of powerful energy stor- 
age which is charged extremely slowly and in a very inefficient way (because 
most of the charging energy is wasted via the energy cascade). For exam- 
ple, if the forcing is concentrated in a thin spherical shell with wavenumber 
lengths between k and k + Ak then the final energy at large scales (assuming 
the energy equipartition) will be k/ Ak times greater than the energy at the 
forcing scales. 

The mechanism of generation of large-scale magnetic fields via nonlinear 
transfer from an energy source at smaller scales may be relevant, in a very 
quahtative sense, to maintenance of large-scale galactic fields by small-scale 
sources provided by supernovae events. To pursue this line further, however, 
one would have to consider a more realistic thin disk geometry, instead of 
a simple homogeneous external field considered in the present paper, which 
could be an interesting future project. 
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